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Abstract. One of the key socioeconomic phenomena to explain is the distribution 
of wealth. Bouchaud and Mezard have proposed an interesting model of economy 
[Bouchaud and Mezard (2000)] based on trade and investments of agents. In the 
mean-field approximation, the model produces a stationary wealth distribution with 
a power-law tail. In this paper we examine characteristic time scales of the model 
and show that for any finite number of agents, the validity of the mean-field result 
is time-limited and the model in fact has no stationary wealth distribution. Further 
analysis suggests that for heterogeneous agents, the limitations are even stronger. We 
conclude with general implications of the presented results. 
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1. Introduction 

Many empirical studies report broad distributions of income and wealth of individuals 
and these distributions are often claimed to have power-law tails with exponents around 
two for most countries [H El El HI E] . The first models attempting to explain the observed 
properties appeared over fifty years ago [HI El E] • Much more recently, physics- motivated 
kinetic models based on random pairwise exchanges of wealth by agents have attracted 
considerable interest P, [TUl [HI [121 [IS]- An alternative point of view is adopted in the 
wealth redistribution model (WRM) where agents continuously exchange wealth in the 
presence of noise [m [151 [IS]- There are also several specific effects which can lead 
to broad wealth distributions [171 [IHl [IS]- (For reviews of power laws in wealth and 
income distributions see [201 [211 [22] , while for general reviews of power laws in science 



In this paper we analyze the WRM with two complementary goals in mind. Firstly 
we investigate the simplest case when exchanges of all agents are identical, focusing 
on the validity of the mean-field approximation which is the standard tool to solve the 
model and derive the stationary wealth distribution. In particular, we show that for any 
finite number of agents there is no such stationary distribution (other finite-size effects 
are discussed for a similar model in [IB]). Secondly we investigate the model's behaviour 
when the network of agent exchanges is heterogeneous. Previous attempts to investigate 
the influence of network topology on the model [HI I25l [23 [21] were all based on the 
mean-field approximation. We show that this is questionable because heterogeneity of 
the exchange network strongly limits the validity of results obtained using the mean-field 
approximation. 

2. Model and its mean field solution 

Adopting the notation used in [14j, we study a simple model of an economy which is 
composed of agents with wealth Vi {i = 1, . . . , N). The agents are allowed to mutually 
exchange their wealth (representing trade) and they are also subject to multiplicative 
noise (representing speculative investments). The time evolution of agents' wealth is 
given by the system of stochastic differential equations (SDEs) 



where a > controls the noise strength. The coefficient Jij quantifies the proportion 
of the current wealth Vj{t) that agent j spends on the production of agent i per 
unit time. We assume the Ito convention for SDEs and dWi{t) is standard white 
noise [221130]. Hence, denoting averages over reahsations by (■), we have {dWi(t)) = 0, 
{dWi{t) dWj{t)) = 6ijdt, and {vi{t) dWi{t)) = 0. By summing dvi{t) over all agents 
one can see that the average wealth VAit) '■= jf^i=i'^i(t) is not influenced by wealth 
exchanges and obeys the SDK dfyi(t) = Wj(t) dPVj(t). Therefore (df^(t)) = 

and {vA{t)) is constant. For simphcity we assume Vi{0) = 1 (i = 1,...,A^) and thus 



see [231l24j.) 




(1) 
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{vi{t)) = 1 and {vA{t)) = 1. (The influence of tlie initial conditions is discussed in 
Section El) 

Tlie system beliaviour is strongly influenced by the exchange coefficients Jij. The 
simplest choice is Jij = J/ {N — 1) where all exchanges are equally intensive — we say 
that the exchange network is homogeneous. By rescaling the time we can set J = 1 
which means that during unit time agents exchange all their wealth. Consequently, ([T]) 
simplifies to 

dvi (t) = {di {t)-Vi{t))dt + V2av, (t) dW, (t) (2) 

where Vi{t) := is the average wealth of all agents but agent i. In the limit 

N ^ oo, fluctuations of Vi{t) are negligible and one can replace Vi{t) — > {vi{t)) = 1 
as in [1A\. Agents then effectively interact only with the "mean fleld" and their wealth 
levels are independent. Using the Fokker-Planck equation for the wealth distribution 
f{vi,t), the stationary solution /(fj) can be found in the form 

A := 1 + (3) 

For fj ^ A — 1, f{vi) decays approximately as a power-law with exponent 2 + 
while the cumulative distribution has exponent 1 + 1/cr^. When Vi is well described by 
([3]), we say that the system is in the power-law regime. 

The empirical studies mentioned above report power-law exponents around 2, 
indicating that in this model, a ~ 1 is needed to obtain realistic power-law behaviour 
of the wealth distribution. In our analytical calculations we assume a < 1; strong noise 
(cr > 1) is discussed separately at the end of the following section. 



f{v^ 



V(\) 



exp 



3. Complete exchange network for a finite 

To examine when the power-law regime is realised, we flrst investigate the time needed 
to reach the mean- fleld solution (j3]). Such relaxation times were studied very recently 
in kinetic models of wealth distribution [28] . 

Given the homogeneous initial conditions Vi{0) = 1 {i = 1, . . . ,N), the exchange 
terms proportional to Vi — Vi are zero at t = and can be neglected for small times. 
Hence when t is small, each Vi{t) evolves independently due to multiplicative noise, Vi{t) 
is lognormally distributed, and its variance is var[fi](i:) = exppcr^t] — 1 = 2a'^t + O(t^): 
we say that the system is in the free regime. From the known variance — cr^) of the 
mean-fleld solution ([3]), we can estimate the transition time ti between the free regime 
and the power-law regime as 



When t S> ti, the system has been given enough time to reach the power-law regime. 

We now recall the average wealth fyi(t). While (df^(t)) = 0, one can see that 
{dv\{t)) is always positive. Hence the variance var[fA(t)] grows without limit, in 
contradiction with the variance of which is flnite for cr < 1. To resolve this 
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disagreement we have to accept that /(fj) as given by ([S]) is not a stationary solution. 
But what comes after the power-law regime? Since the Fokker-Planck equation for 
the joint probabihty distribution f{vi, . . . , vj\f) cannot be solved analytically, we answer 
this question by investigating the average quantities {vf{t)) and {vi{t)vj{t)) {i ^ j); 
now we are considering a < 1 and hence both are well defined. Due to the assumed 
homogeneous network of interactions and the chosen initial conditions, all averages 
{vf{t)) are identical and the same holds for the cross-terms {vi(t)vj{t)); effectively we 
are left with only two variables. From the Ito lemma it follows that d{vf) = (2f j+df j) df , 
and d{viVj) = Vi dvj + Vj dvi + dvi dvj. After substitution of ([2]) and averaging over all 
possible realisations, we obtain the exact set of equations 



dt 



{v.it)v,it))-il-a'){vm 



d{v,{t)v,{t)) _ 2 



(5) 



dt N- 

Since we set fi(0) = 1 (z = 1,...,A^), {vi{t)) = 1 and the initial conditions are 
{vf{0)) = 1 and (fj(O)fj(O)) = 1; for the general case see Section I^?T1 Independently 
of the initial conditions, for cr > 0, has only the trivial stationary solution 
(vf) = (viVj) = 0. This confirms that for a finite A^, there is no stationary distribution 

By solving ([5]) one obtains the variance var[i;j](t) = {vf{t)) — {vi(t))^ as a function 
of time and as a by-product also the correlation between agents i and j 

^/var[vi{t)]vaT[vj{t)] 

Since the resulting expressions are rather complicated, here we discuss only their limiting 
cases. Small time expansions can be found in the form 

vaT[vi]{t) = 2(xH + 0{t^), (7) 

C^J{t) =^t + 0(t^). (8) 

As can be seen, ([71) agrees with our previous reasoning about the log-normal nature of 
f{vi) in the free regime, while ([H]) confirms that in the limit A^ —>■ oo, wealth correlations 
vanish. 

In the limit of large time we obtain 

hm Q,{t) = l-a' + 1 + 0{1/N'). (9) 

Thus, as t increases, the system passes to the synchronized regime where the wealth of 
agents is strongly correlated. One can estimate the transition time by comparing the 
initial linear growth of Cij with its stationary value, leading to 

t2 = {l-(r^)N + 0{l). (10) 

An alternative estimate can be obtained from var[fj](t). Apart from a constant, it 

contains only terms proportional to exp[Ai^2^] where 

_-a^- N{1 - a^) ± - ct2)2 + 2Na^{3 - a^) - a^{A - a^) 

~ 2(Ar- 1) ■ 
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Figure 1. Time evolution of va.r[vi]{t) and Cij{t) for 10 agents (a) and for 10^ agents 
(b). Analytical results following from ([5|) are shown as lines, numerical results obtained 
by averaging over 10^ realisations are shown as symbols, — 0.5. Vertical dotted lines 
indicate the transition times ti, t2, and <3, left to right, respectively. 



Since Ai < 0, A2 > 0, and for a < 1 is |Ai| S> IA2I, the terms proportional to exp[Ait] 
cause the initial saturation of var[t>j](t) but the terms proportional to exp[A2t] eventually 
take over and cause the divergence of var[fj](t). The corresponding transition time can 
be roughly estimated by solving X2t = 1, yielding 

ts = ^^N + 0{l). (11) 

Both ^2 and ^3 describe the transition between the power-law and synchronized regimes: 
the former focuses on the growth of correlations, the latter on the growth of variances. 

To verify the presented analytical results we investigated the model numerically. For 
numerical solutions of stochastic differential equations we used Milstein's method [29| 
|3T] : random numbers were generated using the standard GSL library and the Mersenne 
Twister generator [22], and the time increment was 10~^ in all simulations. In the 
used discretisation scheme, there is a non-zero probability that the wealth Vi{t) becomes 
negative [S3]- However, thanks to the typical value of Vi{t) and the small time step, in 
the presented numerical simulations this was not an issue. As can be seen in Figure [1^, 
our analytical results agree with numerical simulations of the system. Due to the small 
number of agents, transition times ^2,3 are small and the system goes directly from the 
free regime to the synchronized regime. In Figure Wp the number of agents is large 
and the system behaviour is more complex. In the initial period the variance is small 
and correlations are negligible, while in the period t G [2; 300] the variance is almost 
constant and correlations are still small — the system is in the power-law regime (due to 
large computational complexity, no numerical results are shown here). Eventually, for 
t > 10^, the synchronized regime is established. The transition times given by (jl]), ffTOj) . 
and ffTTl) are shown as vertical dotted lines and agree well with the described changes 
of the system behaviour. 

We should sound here a note of caution about the interpretation of the averages 
{vf{t)) and {vi{t)vj{t)) and the wealth distribution f{vi,t). All these quantities are 
ensemble-based: if many copies of the system evolve independently for time t, by 
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Figure 2. Probability density of the average wealth va at various times. Parameter 
values are iV = 10 and — 0.5, probabilities were obtained from 10^ independent 
realisations of the model. 



examining the final wealths of agent i one can estimate both the distribution f{vi,t) 
and the averages. By contrast, when one speaks about an empirical wealth distribution, 
that is based on the wealth of all agents in one realisation only, it is population-based. 
However, when the number of realisations and the number of agents are large and 
the wealth correlations are small, ensemble- and population-based quantities are alike. 
Such behaviour was observable also in the numerical simulations presented above. In 
the free and power-law regimes, the variance of wealth in each realisation was similar to 
var[t'j(t)] (at various times, differences were less than 20% for = 10 and less than 1% 
for = 10 000) and its relative fluctuations between realisations were approximately 
50% for N = 10 and 2% for N = 10 000. As time goes on, fluctuations of the 
population-based variance grow and so does the difference between the ensemble-based 
and population-based variance of wealth. In the synchronized regime, the equivalence 
of the two quantities breaks entirely. 

The nature of the synchronized regime can be better understood by recalling the 
average wealth VA{t) again. As explained above, its evolution is given by a sum of 
multiphcative processes, dvAit) = -^^^iLi^iit) dWi{t). Despite this summation of 
contributions and their variable strengths {{vf{t)) increases with time), according to 
Figure [2] the distribution of va is approximately lognormal and, in agreement with our 
expectations, the variance var[fyi] is increasing. In the initial regime, this increase is 
due to growing variances of all agents' wealth. In the power-law regime, variances 
of wealth levels are approximately constant but their growing correlations lead to 
increasing var[fA]. In the synchronized regime, wealth correlations are already saturated 
and the growth of var[f^] is caused by exponentially growing variances of wealths. 
Since correlations are large, ensemble- and population-based quantities are no longer 
equivalent. Finally we remark that since va > 0, {vA(t)) = 1 is fixed, and var[t>^(t)] 
grows without bounds, in the course of time it is increasingly probable that VA(t) is 
much smaller than its expected value {vA(t)) = 1; this can be interpreted as a high 
occurrence of temporal depressions of the economy. 
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When 0" > 1, both var[fj](t) and Cij{t) diverge and must be replaced by different 
quantities. Instead of the variance, one can use the mean absolute deviation {\vi{t) — l\)'^ 
which avoids second moments of the wealth distribution and hence can be used for any 
(J. The Pearson's correlation coefficient can be replaced by a rank correlation coefficient 
(Kendall's r or Spearman's p). All three proposed quantities are hard to handle in 
analytical calculations and with strong noise, numerical simulations of the system are 
extremely time-demanding. While we have obtained no definite results yet, preliminary 
outcomes suggest that in this case too the transition from the power-law regime occurs 
at a time proportional to the number of agents N . 



4. General exchange network 

Now we generalize the exchange network to an arbitrary graph: denoting the set of 
neighbours of agent i by jV,, the number of neighbours by /c,, the average number of 
neighbours by z. We assume that each agent interacts equally with all neighbours and 
per unit time exchanges the whole wealth, hence 

Jij = 1/kj for i e Mj, Jij =0 for i ^ Mj; (12) 

notice that the matrix of exchanges J is asymmetric. Now, ([T]) generalizes to 

dvi = {vi - Vi) dt + V2avi dWi (13) 

where Vi := '^j^. vj/kj. By averaging over realisations we obtain the set of equations 
for the stationary values of the average wealths 

{«.) = E %^ (") 



mi 



which is solved by {vi) ~ ki. Assuming average wealth equal to 1, ( fT^ has the unique 
solution (f j) = ki/ z. This means that the topology of the exchange network is crucial for 
the distribution of wealth among the agents. Consequently, when a is small and hence 
wealth fluctuations are negligible, a power-law distribution of wealth can be purely a 
topological effect of a scale-free degree distribution in the network of agent exchanges. To 
proceed, {vf(t)) and {vi(t)vj(t)) are again the key quantities. They fulflll the equations 

dM)^ 2$: <f^-2(i-o{«?), 

d{ViVj} _ \^ {VjVm) , \- {ViVn) ^, . 

dt ~ ^ k ^ k 

which can be derived similarly to ([5]). We set the initial conditions according to the 
stationary wealths as fi(0) = ki/z and thus {vf{Q)) = kf/z"^ and (fj(O)fj(O)) = kikj/z"^ 
(the general case is studied in Section HTTI) . From (ITHIl follows 



dt 



>0, 



= 0, (16) 

t=o 
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Figure 3. Time evolution of correlations for the ring network of 10 agents, — 0.25. 
Symbols show numerical results for neighbouring agents (circles), agents with the 
distance 2 (squares), and agents with the distance 3 (diamonds), averaged over 10^ 
realisations. Dashed lines have slopes 1, 2, and 3, respectively. 



which means that the growth of var[fj](t) precedes the growth of Cij{t). This gives us 
a way to investigate the small time behaviour of (1151) : assuming {vi{t)vj{t)) constant, 
we obtain {vf{t)) which in turn leads to an enhanced estimate of {vi{t)vj{t)) . For 
neighbouring agents i and j, the results are 

var[t;,(t)] = ^At + 0(t^), (17) 

a,it) =^^t + oit^). (18) 

Moreover, it can be shown that when the shortest path between agents i and j has 
the length L, the leading term of Cij(t) is proportional to t^. These results are 
confirmed by Figure [3] where we investigate a system of ten agents who are placed 
on a ring {i.e., ki = 2, i = 1, . . . , 10). As can be seen, numerical results agree well with 
Cijit) proportional to t^. The system produces a "cascade" of correlations: first only 
neighbouring agents are considerably correlated, then also agents with the distance two, 
distance three, and so forth. 

For ( JT3l) . the mean-field approximation yields the stationary distribution 

where A'^ = 1 + ki/{za^) and the corresponding variance is var[t>j] = k^a'^/{z'^{l — cr^)). 
By comparing this stationary variance with (|T7I) . we obtain the transition time from the 
free regime to the power-law regime as 



which is identical to (jl]). Further, from ( fTSi) we see that the transition time from the 
power-law regime to the synchronized regime is proportional to kikj/{ki + kj) and thus 
for the whole network it can be estimated as 



4 = Oiz) 



(21) 
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which is a generahzation of ffTUl) . We see that for networks with a relatively small average 
degree, the power-law regime appears only for a limited time or not at all. 

We were unable to obtain an equivalent of the transition time ^3 for a general 
network. Considering, for example, a simple star-like structure with one agent in the 
center and the remaining — 1 agents connected only to him, one can see that the 
transition time ^3 is small and does not scale with N. This suggests that similarly to ^2, 
^3 is also of the order 0{z). This contradicts the findings presented in [Tl] (page 541) 
where they report stationary power-law tails for 2; = 4; it is possible that their numerical 
results are influenced by finite-time and finite-size effects. 

4.I. Influence of the initial conditions 

There is still one more transition time to investigate. When the initial conditions fi(0) 
are not set in line with the stationary wealths given by Eq. ( |T^ . a certain time is 
needed to redistribute the excessive wealth levels over the network; we say that the 
system is in the equilibration regime. Since (dWi) = 0, noise terms do not contribute 
to the redistribution. Thus, ffT^ effectively simplifies to dvi = {vi — Vi) dt which leads 
to the exponential convergence of Vi to the stationary value ki/z. By the substitution 
Ui := Vi — ki/z we obtain 



whose time scale is given by the initial terms iii + Ui as 0{1). Thus, the initial wealth 
distribution equilibrates in time 0(1). Since the transition from the free regime occurs 
roughly at the same time, the system passes from the equilibration regime directly to 
the power-law regime. 

5. Conclusion 

We have shown that in the investigated model, agent wealths have no stationary 
distribution and the power-law tailed distribution reported in previous works is only 
transient. In addition, for any finite number of agents, their average wealth va follows 
a multiplicative process with a fixed expected value (va) and an increasing variance 
var[v^]. Hence, as illustrated in Figure [2l the probability P{va < x) approaches 1 for 
any x > 0. We can conclude that the simple economy produced by the model is an 
uneasy one: the longer it evolves, the higher the probability that a given agent i has 
wealth much smaller than any positive fraction of the expected wealth (fj). 

There is also a more general lesson to be learned. In essence, the mean-field 
approximation here anchors the agent wealths to their expected values and thus weakens 
the diffusive nature of the studied stochastic system. Mathematically speaking, the 
system behaviour depends on the order of limits N ^ 00 and t ^ 00: in the former 
case there is a stationary wealth distribution, in the latter case there is none. This is an 
undesired consequence of the mean-field approximation which, as with other stochastic 




(22) 
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models, should be used with great caution. In particular, when using it, one should 
check if the nature of the studied system is not changed. To achieve this, in this paper 
we have used an aggregate quantity (the average wealth) and a quantity obtained using 
the mean- field approximation (the wealth variance). 

On the other hand, in some cases an anchoring term may be appropriate. For 
example, a simple taxation of wealth can be achieved by introducing the term r(l — fj) dt 
to ([2]), where r > represents the tax rate. Then the set of equations for {vf{t)) and 
{vi{t)vj(t)) has a nontrivial stationary solution for a < 1; one can say that the proposed 
taxation stabilizes the system. Notably, systems of coupled stochastic equations with 
multiplicative noise and negative feedback are common in the study of nonequilibrium 
phase transitions in magnetic systems [M] . Our work shows that this negative is crucial 
for mean-field studies of such systems [35] . 

In addition to the presented results, several questions remain open. First, for large 
time t, the analytical form of the wealth distribution f{vi, t) is unknown. Second, for an 
arbitrary network of exchanges, the limiting value of the correlation Cij{t) and also the 
transition time t'^ are of interest. Third, the strong noise case deserves more attention 
and perhaps an attempt for approximate analytical results. Finally, the studied model 
is simplistic, since it combines two ingredients of economy — trade and speculation — in 
a very unrealistic way. Devising a more adequate model remains a future challenge. 



Acknowledgments 

We acknowledge the hospitality of the Comenius University (Bratislava, Slovakia) and 
the Fribourg University (Fribourg, Switzerland). We thank Frantisek Slanina for early 
discussions and Zoltan Kuscsik and Joseph Wakeling for helpful comments. 



References 

[1] Pareto V 1897 Cours d'economie politique (Lausanne: Rouge) 
[2] Piggott J 1984 Economic Record 60 252-265 
[3] Aoyama H et al 2000 Fractals 8 293-300 

[4] Dragulescu A and Yakovenko V M 2001 Ptiysica A 299 213-221 
[5] Sinha S 2006 Physica A 359 555-562 

[6] Champernownc D G 1953 Ttie Economic Journal 63 318-351 

[7] Wold H O A and Whittle P 1957 Econometrica 25 591-595 

[8] Stiglitz J E 1969 Econometrica 37 382-397 

[9] Ispolatov S et al 1998 Eur. Phys. J. B 2 267-276 
[10] Dragulescu A and Yakovenko V M 2000 Eur. Phys. J. B 17 723-729 
[11] Slanina F 2004 Pliys. Rev. E 69 046102 

[12] Patriarca M et al 2006 The ABCD's of statistical many-agent economy models Preprint 

|arXiv:physics/0611245 
[13] Chatterjee A and Chakrabarti B K 2007 Eur. Phys. J. B 60 135-149 
[14] Bouchaud J-P and Mezard M 2000 Physica A 282 536-545 
[15] Solomon S and Richmond P 2001 Physica A 299 188-197 

[16] Di Matteo T et al 2004 in The Physics of Complex Systems (New Advances and Perspectives) Eds. 
Mallamace F and Stanley H E (Amsterdam: lOS Press) 



Breakdown of the mean-field approximation in a wealth distribution model 



11 



Sornette D 1998 Phys. Rev. E 57 4811-4813 

Huang Z-F and Solomon S 2001 Physica A 294 503-513 

Reed W J 2001 Economics Letters 74 15-19 

Quadrini V and Rios-Rull J-V 1997 Models of the distribution of wealth Research Department, 

Federal Reserve Bank of Minneapolis 
Davies J B and Shorrocks A F 2000 The Distribution of Wealth Handbook of Income Distribution 

(Handbooks in Economics) Eds. Atkinson A B and Bourguignon F (Amsterdam: North Holland) 
Yakovenko V M 2007 Statistical Mechanics Approach to Econophysics Springer Encyclopedia of 

Complexity and System Science (Springer Verlag) 
Newman M E J 2005 Contemporary Pfiysics 46 323-351 

Farmer J D and Geanakoplos J 2006 Power laws in economics and elsewhere SFI Technical report 



Souma W et al 2001 Small- World Effects in Wealth Distribution Preprint |cond-mat/0l08482] 
GarlascheUi D and Loffredo M I 2004 Ptiysica A 338 113-118 
GarlascheUi D and Loffredo M I 2008 J. Phys. A: Math. Theor. 41 224018 
Gupta A K 2008 Physica A 387 6819-6824 

Gardiner C W 2004 Handbook of Stochastic Methods, 3rd Edition (Berlin: Springer) 

Van Kampen N G 2007 Stochastic Processes in Physics and Chemistry, 3rd Edition (Amsterdam: 

North Holland) 
Higham D J 2001 SIAM REVIEW 43 525-546 

Galasi M et al 2002 CNU Scientific Library (Reference Manual) (Bristol: Network Theory Ltd) 
Moro E 2004 Phys. Rev. E 70 045102(R) 
Birner T et al 2002 Phys. Rev. E 65 046110 
Mufioz MA et al 2005 Phys. Rev. E 72 056102 



